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ABSTRACT 

We present the Chandra ACIS and ASCA CIS results for a series of four long-term observations 
on DoAr 21, ROXs 21 and ROXs 31; the X-ray brightest T-Tauri stars (TTSs) in the p Ophiuchi 
cloud. In the four observations with a net exposure of ~ 600 ksec, we found six, three and 
two flares from DoAr 21, ROXs 21 and ROXs 31, respectively; hence the flare rate is fairly 
high. The spectra of DoAr 21 are well fitted with a single-temperature plasma model, while 
those of ROXs 21 and ROXs 31 need an additional soft plasma component. Since DoAr 21 is 
younger (~10 5 yr) than ROXs 21 and ROXs 31 (~T0 6 yr), these results may indicate that the soft 
component gradually increases as T-Tauri stars age. The abundances are generally sub-solar and 
vary from element to element. Both high-FIP (first ionization potential) and low-FIP elements 
show enhancement over the mean abundances. An unusual giant flare is detected from ROXs 
31. The peak luminosity and temperature are ~10 33 ergs s _1 and ~10 keV, respectively. The 
temperature reaches its peak value before the flux maximum, and is nearly constant (4-5 keV) 
during the decay phase, indicating successive energy release during the flare. The abundances 
and absorption show dramatic variability from the quiescent to flare phase. 
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1. INTRODUCTION 

Low-mass stars in the pre-main sequence 
(PMS) stage are classified into four types based 
on their infrared (IR) to sub-mm spectral en- 
ergy distributions (SEDs); molecular cloud cores 
through protostars are represented by class and I 
SEDs, while classical and weak-line T-Tauri stars 
(CTTSs and WTTSs) exhibit class II and III SEDs 
(Shu, Adams, & Lizano 1987; Andre & Montmerle 
1994). 

The Einstein satellite first revealed that T- 
Tauri stars (TTSs = CTTSs and WTTSs) emit 
intense and time variable X-rays with occasional 
rapid flares (Feigelson & DeCampli 1981; Mont- 
merle et al. 1983). These properties are consistent 
with a scenario of enhanced solar-type activity, at- 
tributable to magnetic dynamo processes. Succes- 
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sive X-ray satellites such as ROSAT, ASCA, and 
Chandra have detected X-ray emission from many 
PMS stars, even from the youngest class stage 
(Tsuboi et al. 2001), and suggest that X-ray emis- 
sion from low-mass PMS stars is caused by mag- 
netic activity. These results, however, rely mainly 
on short duration observations; neither long-term 
nor detailed behavior of the X-ray spectrum and 
timing has been well studied. This paper hence ad- 
dresses three topics on TTSs; (1) the duty ratio of 
the X-ray flare, (2) the long-term plasma structure 
and variability, and (3) the chemical compositions 
in the quiescent and flare phases. 

(1) Stelzer, Neuhauser, & Hambaryan (2000) 
compiled ROSAT PSPC observations of the 
Taurus- Auriga-Perseus region and found that the 
flare rate may depend on the stellar age, rotation 
period, and/or multiplicity; stars with younger 
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ages, shorter rotation periods, and located in bi- 
nary systems tend to show more frequent flares. 
Hence, the flare rate may help to reveal the phys- 
ical conditions of PMS stars. 

(2) Prcibisch (1997) fitted the ROSAT PSPC 
spectra of ^60 PMS stars with several thermal 
plasma models and found that multi-temperature 
components are needed for most of the spec- 
tra. Tsujimoto et al. (2001) also reported that 
about two-thirds of the Chandra spectra of TTSs 
in the Orion molecular cloud tend to be well 
fitted by multi-temperature plasma models as 
compared with single temperature plasma mod- 
els for younger objects (protostars). The multi- 
temperature structure may also depend on the 
activity phase of each star; higher temperature 
plasma becomes dominant in a flare compared 
with a quiescent phase. Thus long-term monitor- 
ing of both flare and quiescent phases is crucial 
for the study of the origin and structure of the 
multi-temperature plasma. 

(3) ASCA found that the mean chemical com- 
position of X-ray-emitting plasma in young stellar 
objects (YSOs) is generally sub-solar and varies 
from star to star (Kamata et al. 1997; Ozawa 
2000). These results may be misleading, however, 
because the mean-abundance determination is bi- 
ased toward elements with strong emission lines 
in the relevant energy range and largely depends 
on the plasma temperature. Moreover, the abun- 
dances of the stellar corona may be different from 
those in the photosphere, depending on the first 
ionization potential (FIP) of the relevant elements. 
For instance, in a solar corona, elements with low 
FIP are overabundant relative to the solar photo- 
spheric values ("the FIP effect", Feldman 1992). 
In some of the stellar coronae, however, the abun- 
dances of high-FIP elements are larger than the 
others (the Inverse FIP = IFIP effect). The IFIP 
effect appears in X-ray active stars, while the FIP 
effect is generally seen in less active stars (Au- 
dard, Giidel, & Mewe 2001; Brinkmanet al. 2001; 
Gudel et al. 2001a, b,c; Huenemoerder, Canizares, 
& Schulz 2001; Kastner et al. 2001). The depen- 
dence of abundance on the X-ray activity is also 
found within individual stars; the enhancements in 
abundance are seen from quiescent to flare phases, 
or even within a single flare phase (e.g. HR1099, 
Audard et al. 2001; V773 Tau, Tsuboi et al. 1998). 
Due to the limited statistics and samples, however, 



it is unclear whether or not abundance variations 
between quiescent and flare phases are common, 
and whether or not FIP/IFIP effects are present 
in low-mass YSOs. 

To investigate the above-mentioned subjects, 
we selected the three X-ray brightest TTSs; DoAr 
21 (ROXs 8), ROXs 21 and ROXs 31 in the p Ophi- 
uchi cloud (p Oph) located at a distance of 145 pc 
(de Zeeuw et al. 1999), because all three have pub- 
licly available long-term observations with Chan- 
dra and ASCA. These TTSs were first recognized 
as bright X-ray sources with Einstein (Montmerle 
et al. 1983), and were identified with optical K- 
M stars having weak Ha emission (Bouvier & Ap- 
penzeller 1992), hence classified as WTTS (class 
III). Simon et al. (1995) observed the lunar oc- 
cultation of these WTTSs at infrared wavelengths 
and found that ROXs 21 and ROXs 31 are binary 
systems with separation angles of 0'.'3 and 0'.'48, 
respectively. The rotational period of ROXs 21 
is reported to be 1.39 days (Shevchenko & Herbst 
1998). Since no continuum emission is found at 1.3 
mm (Andre & Montmerle 1994), the circumstellar 
envelope of these TTSs is thought to have already 
disappeared, while DoAr 21 may still have an ac- 
cretion disk, as suggested by the detection of near- 
infrared (NIR) polarized emission (Ageorges et al. 
1997). Using the theoretical evolutionary tracks 
in the H-R diagram, the age of DoAr 21 is esti- 
mated to be ~10 5 yr, which is younger than that 
of ROXs 21 and ROXs 31 (-10 6 yr) (Niirnberger 
et al. 1998). Finally, DoAr 21 and ROXs 31 show 
variable radio emission, which is probably due to 
the gyro-synchrotron emission induced by the sur- 
face magnetic field (Stine et al. 1988). 

Early results with the Chandra and ASCA ob- 
servations (Koyama et al. 1994; Kamata et al. 
1997; Tsuboi 1999; Tsuboi et al. 2000; Skin- 
ner 2000; Imanishi et al. 2001) revealed that 
these particular TTSs have luminosities of >10 30 
ergs s _1 and relatively low absorption (Ah < 10 22 
cm -2 ). Thus they can provide enough counts to 
measure their temporal and spectral properties in 
great detail. 

2. OBSERVATIONS AND DATA RE- 
DUCTION 

The Chandra (Weisskopf, O'dcll, & van Spey- 
broeck 1996) observation (hereafter, obs CI) of 
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these TTSs was made with five front-illuminated 
X-ray CCDs, ACIS-IO, II, 12, 13, and S2, in the 
0.2-10 keV band. The level 2 data were retrieved 
from the Chandra X-ray Center (CXC) archive. X- 
ray events are selected using the ASCA grades 0, 

2, 3, 4, and 6. An effective exposure time of «100 
ks is obtained (Table 1). Another long exposure 
observation of the p Oph core A, which contains 
DoAr 21 in the ACIS-I3, is also in the archive (obs 
C2), but this data suffer from the photon pile-up 
and therefore were not used. 

Three observations (obs Al, A2, and A3) were 
carried out with the two Gas Imaging Spectrom- 
eters (GISs) (Ohashi et al. 1996) and the two 
Solid-state Imaging Spectrometers (SISs) (Burke 
et al. 1991) onboard ASCA (Tanaka, Inoue, & 
Holt 1994), at the foci of the X-ray telescopes 
(XRTs) (Scrlemitsos et al. 1995) sensitive to pho- 
tons in 0.4-10 keV. However, the TTS sources 
were located outside of, or at the edge of the SIS's 
field of view in the majority of the observations; 
therefore we do not use the SIS data. We retrieved 
the ASCA unscreened data from the HEASARC 
Online Service. Data taken at a geomagnetic cut- 
off rigidity lower than 4 GV, at an elevation angle 
less than 5° from the Earth, and during passage 
through the South Atlantic Anomaly were filtered 
out. Particle events were also removed using the 
rise-time discrimination method. The total avail- 
able exposure times were «38 ks (obs Al), 93 ks 
(A2), and 75 ks (A3) (Table 1). 

3. ANALYSIS AND RESULTS 

3.1. Extraction of Source and Background 
Photons 

The positions of NIR counterparts (Barsony et 
al. 1997) were used to fine-tune the Chandra co- 
ordinates. We then extracted the source photons 
using an ellipse with major and minor axes of 60'.'0 
and 35'.'0, respectively, for DoAr 21, and from re- 
spective circles of 11'.' 7 and 12'.'0 radius for ROXs 
21 and ROXs31, with the regions centered at the 
NIR positions (Barsony et al. 1997). The rela- 
tively large source areas are due to large off-axis 
point spread functions (PSFs); DoAr 21 lies in 
ACIS-S2 ~20' away from the optical axis, while 
ROXs 21 and ROXs 31 are at the edge of ACIS- 
II with off-axis angles of ~7-10'. The large off- 
axis angle (and hence large PSF) for the bright- 



est source DoAr21 mitigates photon pile-up, which 
in fact is a problem for the obs C2 case. The 
background region for DoAr 21 was taken from a 
source-free 19 arcmin 2 region on ACIS-S2, while 
that for ROXs 21 and ROXs 31 was taken from 
a 63 arcmin 2 region on the ACIS-I array. For the 
ASCA GIS data (obs A1-A3) the source photons 
were extracted from 3' radii circles. The back- 
ground regions were also taken using 3' radii cir- 
cles from nearby source-free regions. 

3.2. Time Variability 

Figure 1 and 2 show the X-ray light curves (no 
background subtraction) of (a) DoAr 21, (b) ROXs 
21, and (c) ROXs 31 obtained with the Chandra 
ACIS (obs CI) and the ASCA GIS (obs A1-A3), 
respectively. The photon counts are normalized 
by the effective area at 1 keV, for comparison with 
the different instruments and observation epochs. 
The time bins of each data point are 10 3 s and 
10 4 s for ACIS (Figure 1) and GIS (Figure 2), 
respectively. The energy band of the ACIS and 
GIS data for DoAr 21 arc taken to be 0.5-9.0 keV. 
However, the GIS light curves for ROXs 21 and 
ROXs 31 are limited in the soft energy band of 0.5- 
1.5 keV in order to avoid possible contamination 
from a nearby variable hard X-ray source, YLW 
15A(Tsuboi 1999; Tsuboi et al. 2000; Imanishi et 
al. 2001). To verify the energy band, we extracted 
the ACIS-I spectrum of YLW 15A and fit it with 
an absorbed thin thermal plasma model. We find 
a best-fitted absorption of Ah = 4.4xl0 22 cm~ 2 
(Imanishi et al. 2001), which reduces the X-ray 
emission below 1.5 keV to less than 3% of the 1.5- 
9 keV flux. Hence possible contamination from 
YLW 15 A to the GIS source circles of ROXs 21 
and ROXs 31 is insignificant in the 0.5-1.5 keV 
band. No data is shown in the ROXs 21 light curve 
(Figure 2b) from MJD = 50511.55-50512.20 due 
to significant contamination from an exceptionally 
giant flare of ROXs 31 (see Figure 2c). 

For a unified study, we conventionally define 
a flare using the following criterion: at least two 
consecutive time bins must have larger fluxes than 
the quiescent bin by > 5cr-level, where the quies- 
cent bin is the minimum time bin in each obser- 
vation. Under this criterion, we detect six, three, 
and two flares from DoAr 21, ROXs 21 and ROXs 
31, respectively. These are labelled as "F", "Fl", 
"F2", and "F3" with arrows in Figure 1 and 2, 
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while the quiescent phases are indicated by "Q". 
An extremely large flare with the maximum flux of 
~100 times higher than the quiescent level is de- 
tected from ROXs 3f (Figure 2c) and is separately 
treated in §3.4. 

3.3. X-ray Spectra 

The background-subtracted ACIS spectra of 
each source using all the data of obs CI are shown 
in Figure 3. Several emission lines from highly ion- 
ized elements such as O, Ne, Na, Mg, Si, S, Ar, Ca, 
and Fe are seen, hence the spectra are characteris- 
tic of a thin thermal plasma. We therefore fit the 
spectra with a thin thermal plasma model allow- 
ing the abundance of each element with prominent 
lines (indicated by arrows in Figure 3) to be free. 
For the study of time variation of the spectra, we 
separately fit the flare and quiescent spectra as 
shown in Figure 1. A single-temperature (1-T) 
thin-thermal plasma model (Mewe et al. 1985) is 
not acceptable for all the spectra except DoAr 21 
and the flare phase of ROXs 31. We therefore fit 
two-temperature (2-T) plasma models for ROXs 
21 and ROXs 31 and find acceptable fits. Table 2 
shows the best-fit parameters. All the elemental 
abundances are below the solar photosphcric val- 
ues (Anders & Grevesse 1989) except for Ne in 
DoAr 21 and Na in ROXs 21. The best-fit abun- 
dances are displayed in Figure 4 with the order of 
the FIP. 

We investigate the Na overabundance seen in 
ROXs 21 in further detail, because Na is a rela- 
tively rare clement compared with the major 4N 
nuclei and the Na Lyman-a line lies near the 
Lyman-/? line of the more abundant element Ne. 
Figure 5a shows the best-fit result when the Na 
abundance is correctively varied with the "other" 
elements. A significant data excess is found at 
«1.2 keV, near the Lyman- (3 line of Ne (1.211 
keV). The intensity ratio of Lyman-/? to -a lines 
of Ne should increase with increasing plasma tem- 
perature, but stay almost constant at « 7 for 
plasma temperatures higher than 1 keV (Mewe 
et al. 1985). Therefore 1/7 of the Lyman-a flux 
should be an upper-limit to the Lyman-/? line flux. 
Since the flux of the Ne Lyman-a line is fairly 
strong, the Lyman-/? line flux is well constrained 
as seen in Fig 5a. Thus the excess flux should orig- 
inate from another element; the Lyman-a line of 
Na (1.236 keV) is the best candidate. Allowing the 



Na abundance to be a free parameter improved the 
fit as demonstrated in Figure 5b. The reduction 
of x 2 (d.o.f.) from 114 (97) to 108 (96) is signif- 
icant with ~98 % confidence level for the F-test 
(Bevington & Robinson 1992). 

Using the GIS spectra of DoAr 21, we investi- 
gate the long-term (^yr) behavior of the spectral 
parameters. Since the energy resolution does not 
allow us to separate the emission lines from the key 
elements, we correctively varied the abundances. 
The best-fit luminosity shows long-term variabil- 
ity in the range (2 - 8) x 10 31 erg s -1 . The best- 
fit temperatures and abundances also vary from 
2.5 to 4 keV and 0.15 to 0.5 solar, respectively; 
both are possibly correlated to the X-ray luminos- 
ity as are shown in Figure 6a and 6b. The Nu 
value is consistent with being constant, although 
weak variability within a factor w 2 cannot be 
rejected. We note that although the 1-T model 
in the medium resolution X-ray spectra may lead 
to an error in abundance by a factor ~2 (Reale, 
Peres, & Orlando 2001), such an effect is signifi- 
cantly reduced with 2-T or multi-T models. Since 
we used 2-T models for 2 out of the 3 sources, the 
artificial abundance uncertainty should be better 
constrained. 

3.4. Giant Flare from ROXs 31 

For the study of the spectral evolution during 
the flare, we make and fit time-sliced GIS spectra 
of the giant flare from ROXs 31 in obs A2 (Fig- 
ure 2c). The time intervals are shown in Figure 
7 (left). The spectra show strong Ka emission of 
He-like Fe at «6.7 keV, hence we let the Fe abun- 
dance vary. The flare spectra are well fitted with a 
1-T model. The best-fit parameters for each time 
interval are shown in Figure 7 (right). As shown 
in the figure, this flare has an unusual time pro- 
file. The light curve increases slowly over ~15 ks 
(see also Skinner 2000). The temperature reaches 
its peak value before the flux maximum, and stays 
nearly constant during the decay phase. 

The absorption (Ah) and abundances dramat- 
ically change from quiescent to flare; in the 
flare is ~7.1 xlO 21 cm~ 2 , ^2 times smaller than 
that of the quiescent phase (obs CI) of 1.7 xlO 22 
cm~ 2 , and the abundances increase by a factor of 
3-30 from quiescent (Z = 0.01-0.1 Z Q ) to the gi- 
ant flare (0.26-0.29). These parameters, however, 
do not change during the flare. 
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4. DISCUSSION 

4.1. Flare Rate 

We have detected six (DoAr 21), three (ROXs 
21) and two (ROXs 31) flares under our flare cri- 
terion (see §3.2). The total exposure time of the 
four observations is —7 days, hence the flare rate is 
one per 1.2, 2.3, and 3.5 days for DoAr 21, ROXs 
21 and ROXs 31, respectively. We also confirmed 
the high flare rate of DoAr 21 in obs C2, although 
the data suffer from photon pile-up. Two possibil- 
ities could account for the higher rate in DoAr 21: 
differences of the energy band (DoAr 21: 0.5-9.0 
keV, versus ROXs 21 and ROXs 31: 0.5-1.5 keV 
in obs A1-A3) and the mean count rate (DoAr 21: 
-0.4 counts s" 1 , versus ROXs 21 and ROXs 31: 
0.04-0.09 counts s" 1 in obs CI). In fact, limiting 
the GIS data of DoAr 21 to 0.5-1.5 keV reduces 
the number of flares (under our criterion) to two 
(F2 in obs A2 and F in obs A3), similar in num- 
ber to those of ROXs 21 and ROXs 31. Also, the 
higher count rate of DoAr 21 results in a higher 
sensitivity to smaller amplitude flares that can be 
detected under our criterion (§3.2). Indeed, both 
flares in Figure la have smaller amplitudes than 
those in Figure lb and c. Conversely, at similar 
sensitivities, the flare rate of ROXs 21 and ROXs 
31 may be comparable to that of DoAr 21, in spite 
of different ages (DoAr 21, -10 5 yr; ROXs 21 and 
ROXs 31, -10 6 yr) and different structure (ROXs 
21 and ROXs 31 are binaries, while DoAr 21 may 
be a single star). 

The typical flare rate of X-ray sources in the 
Taurus- Auriga-Perseus region (Stelzer et al. 2000) 
is 1/(4-5) days (assuming a typical decay time 
scale of 1 hour) . Therefore we predict significantly 
higher flare rates than that reported previously. 
The higher duty ratio may be primarily due to 
the extended sensitivity in the hard X-ray band 
(>1.5 keV) and their brightness, because the flare 
activity (flux increase) is clearer in the harder X- 
ray band and/or for brighter sources as we have 
already demonstrated for DoAr 21 in the previous 
paragraph and in Figure 6a. 

4.2. Interpretation of the Temperature 
Structure 

The high-quality spectra of the long Chandra 
ACIS exposure reveal that the 2-T models gives a 



better fit for overall spectra than the 1-T models 
for ROXs 2f and ROXs 31. This supports pre- 
vious 2-T model fits for some fraction of other 
bright TTSs (Carkner et al. 1996; Preibisch 1997; 
Costa et al. 2000; Ozawa 2000; Tsujimoto et 
al. 2001). DoAr 21, on the other hand, displays 
a simple 1-T spectrum with kT —3 keV; hence 
it has no additional soft component. Since the 
age of DoAr 21 (~ 10 5 yr) is younger than that of 
ROXs 21 and ROXs 31 (-10 6 yr) (Niirnberger et 
al. 1998), coupled with the result of Tsujimoto 
et al. (2001) that 2-T spectra are found more of- 
ten in older TTSs than in younger protostars, we 
speculate that the soft component, probably a rel- 
atively steady corona, is generated gradually as 
the system increases in age, finally reaching solar- 
like corona. In this scenario, the hard component 
would be the sum of unresolved flares. 

4.3. Chemical Composition of the X-ray- 
emitting Plasma 

The coronal abundances of the TTSs are sub- 
solar, consistent with the previous results. From 
Figure 4, we see that both high-FIP (Ne and Ar) 
and low-FIP (Na, Mg and Ca) elements show 
higher abundances than the other elements (the 
IFIP and FIP effects). 

For the abundances in solar corona, the FIP- 
effect appears in elements with FIPs below 10 eV; 
these elements are collisionally ionized in the pho- 
tosphere at 6000-7000 K temperature, and would 
be preferentially transferred to the upper coronal 
region by electric fields (Feldman 1992). DoAr 
21 and ROXs 21 are a K0 star and a binary of 
K4 with M2.5, respectively, and hence have pho- 
tospheric temperatures of 4000-5000 K, 0.6-0.8 
times that of the solar photosphere (Niirnberger 
et al. 1998). Therefore the FlP-effect energy- 
limit of 10 eV should be shifted to 6-8 eV, which 
is near Mg (FIP = 7.8 eV) and Ca (6.1 eV), but 
well above Na (5.1 keV). Hence, the abundance 
enhancements of Mg, Ca and Na provide indepen- 
dent evidence supporting the universality of the 
FlP-scenario proposed by Feldman (1992). 

The X-ray luminosities of our TTS sample are 
a few x 10 30-31 ergs s" 1 , which is significantly 
higher than typical FIP sources of < 10 29 ergs s _1 
(such as the sun, n 1 UMa, and x 1 Ori; Giidel et al. 
2001c). The high luminosity may also be respon- 
sible for the IFIP effect we observed. A related 
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issue would be the dramatic abundance increase 
in the giant flare of ROXs 31. 

A scenario for the IFIP effect was proposed by 
Gudcl et al. (2001c). A "downward-pointing" 
electric field is produced by the high-energy non- 
thermal electrons during X-ray active phase, and 
suppresses the transfer of the ionized low-FIP ele- 
ments into the corona. The electrons also loose 
their kinetic energy at the photosphere by the 
collisional process and then heat up the cool ele- 
ments, making the abundance increase during the 
flare. The radio gyro-synchrotron emission from 
DoAr 21 and ROXs 31 further supports this idea, 
because this is the direct evidence for the presence 
of the high-energy electrons. Since the abundance 
increase should be correlated with the non-thermal 
radio flux variations, simultaneous radio and X- 
ray observations would help to constrain physical 
models. 

Another possibility is photo-ionization by coro- 
nal X-rays. This was proposed to explain a solar 
flare in which Ne is enhanced relative to O (Shemi 
1991; Schmelz 1993). The IFIP effect is seen in 
the X-ray brightest stars with Lx = 10 29 -10 32 
ergs s" 1 (HR 1099, Brinkman et al. 2001; AB Dor, 
Giidel et al. 2001a; YY Gem, Gudcl et al. 2001b; 
II Pegasi, Hucnemoerder et al. 2001; and TW Hy- 
drae, Kastner et al. 2001), in which large X-ray 
photo-ionization may occur. The abundances in- 
crease in the X-ray flare may also be explained by 
this scenario. 

4.4. Interpretation of the Giant Flare of 
ROXs 31 

The unusual giant flare from ROXs 31 shows 
1-T plasma emission, in contrast to the quiescent 
phase. The flux of the soft component in the qui- 
escent phase is 10-100 times lower than that of 
the flare flux. Hence it is highly possible that the 
soft component does exist during the flare but is 
hidden by the more luminous harder component. 

The peak luminosity of ^10 33 ergs s _1 and tem- 
perature of ~10 keV are comparable to the gi- 
ant flare of another TTS binary system V773 Tau 
(Tsuboi et al. 1998). However, the time pro- 
file of the plasma temperature is different from 
that of V773 Tau; the temperature arrives at the 
peak value preceding the flare peak, and, in the 
decay phase, stays nearly constant at 4-5 keV. 



These imply the existence of successive energy re- 
lease due to the magnetic reconnection and that 
the generation and heating of the plasma occur 
even in the decay phase, as suggested by Tsuboi 
et al. (1998); Schmitt & Favata (1999); Favata 
& Schmitt (1999). Such a phenomenon is also 
seen in long duration event (LDE) flares of the 
Sun (Isobc et al. 2001). 

The abundance and absorption during the gi- 
ant flare are systematically larger and smaller than 
those in the quiescent phase (Table 2). This could 
be due to evaporation of photosphere and envi- 
ronments by high-energy electrons, which may be 
produced just after magnetic reconnection. Evap- 
oration by the flare X-rays may not be important, 
because the abundances and absorption do not 
change much during the flare. Thus, in the case 
of the IFIP effect in the giant flare of ROXs 31, 
the high-energy electron origin is more likely than 
X-ray irradiation scenario. 

5. SUMMARY 

We analyzed the Chandra and ASCA data of 
the three X-ray-bright TTSs, DoAr 21, ROXs 21, 
and ROXs 31. Based on the timing and spectral 
features, we derived several results as follows; 

1. Six flares in the 0.5-9 keV band are detected 
from DoAr 21 during the total duration of 
^7 days. Also, three and two flares are 
detected from ROXs 21 and ROXs 31, re- 
spectively, with Chandra in 0.5-9.0 keV and 
ASCA in 0.5-1.5 keV. Converting to the 0.5- 
9 keV band sensitivity, the flare numbers be- 
come nearly the same as DoAr 21. Thus the 
flare rate of these TTSs may be extremely 
high (about one per one day). 

2. The spectra of DoAr 21 are well fitted with 
a 1-T model, while those of ROXs 21 and 
ROXs 31 need additional soft components. 
This can be interpreted as a gradual genera- 
tion of the soft component with the increas- 
ing age. 

3. We found relative abundance enhancements 
in both the low-FIP (Na, Mg and Ca) and 
high-FIP (Ne and Ar) elements. 

4. An unusual giant flare is detected from 
ROXs 31. The temperature remains high 
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(4-5 keV) even during the decay phase, sug- 
gesting a successive energy release. A large 
abundance and absorption jump from quies- 
cent to flare is found, which may suggest a 
large contribution to gas evaporation from 
the photosphere and environments by non- 
thermal electrons. 
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Fig. 1 — Light curves of (a) DoAr 21, (b) ROXs 21, and (c) ROXs 31 in 0.5-9.0 keV obtained with Chandra 
ACIS (obs CI). The vertical axis is the X-ray count rate normalized in the effective area of 100 cm 2 at 1 
keV. The lower arrows indicate the times of the flare and quiescent phases (see Table 2). 



Al A2 A3 Al A2 A3 Al A2 A3 




MJD (day) MJD (day) MJD (day) 



Fig. 2. - Light curves of (a) DoAr 21 in 0.5-9.0 keV, and of (b) ROXs 21 and (c) ROXs 31 in 0.5-1.5 keV, 
obtained with ASCA GIS (obs A1-A3). The vertical axis is the X-ray count rate normalized in the effective 
area of 100 cm 2 at 1 keV. The lower arrows indicate the times of the flare and quiescent phases. 
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Fig. 3. — Time-averaged X-ray spectra of (a) DoAr 21, (b) ROXs 21, and (c) ROXs 31 obtained with 
Chandra ACIS-I (obs CI). The prominent spectral features are indicated by the arrows and the corresponding 
elements. 
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(a) DoAr 21 



(b) ROXs 21 



(e) ROXs 31 



Fig. 4.— Elemental abundances of (a) DoAr 21, (b) ROXs 21, and (c) ROXs 31. Errors indicate 90% 
confidence limits. The horizontal axis shows corresponding elements with their FIP values (eV). Gray 
squares, black-open, and black-filled circles represent abundances in the quiescent, first and second flare 
phases, respectively. 




1.2 

Energy (keV) 



1.2 

Energy (keV) 



Fig. 5. — Close-up view of the fitting results of ROXs 21 in phase Q, in which the Na abundance (a) is varied 
correctively with the "other" elements and (b) is treated as a free parameter. The positions of K-shell lines 
from (a) Ne and (b) Na are shown by arrows. The prominent feature near the Lyman-a line of Na is also 
indicated by dashed squares. 
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Fig. 6. — Plot of the X-ray luminosity and the plasma temperature (a) and the mean abundance (b) of DoAr 
21. Circles and triangles represent the data of ACIS (obs CI) and of GIS (obs A1-A3). Errors indicate the 
90% confidence limits. 




Fig. 7. — (left) A magnified view of the light curve of ROXs 31 during the giant flare in 0.5-9.0 keV 
(Figure 2c). (right) Time profiles of the best-fit temperature, absorption column, abundances of Fe and 
other elements, and X-ray luminosity of ROXs 31 during the giant flare in obs A2. Errors indicate the 90% 
confidence limits. 
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Table 1 
Observation Log 



Obs.ID 


Detector 


Sequence ID 


Date 


Exposure 
(ks) 


Al 


ASCA/GIS 


20015010 


1993 Aug 20-20 


37.7 


A2 


ASCA/GIS 


25020000 


1997 Mar 2-4 


93.1 


A3 


ASCA/GIS 


96003000 


1998 Aug 13-15 


74.7 


CI 


Chandra/ ACIS 


200060 


2000 Apr 13-14 


100.6 
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Table 2 

Best-Fit Parameters in obs Cl a 



Source 




DoAr 21 






ROXs 21 


ROXs 31 


Phasc b 


o° 

t 


Fl 


F2° 




Q 


F 


Q 


F 


kTi (keV) 


2.5(2.4-2.6) 


2.9(2.7-3.1) 


2.9(2.8-3.1) 





69(0.65-0.74) 


0.86(0.78-4.5) 


0.23(0.15-0.31) 


0.27(0.21-0.41) 


kT 2 (kcV) 








1 


8(1.6-2.1) 


2.6(2.3-3.1) 


1.5(1.2-1.7) 


2.0(1.8-2.5) 


AT /in22 ,. vvl — 2\ 

JVh t-LU cm J 


i n^n (17 i i 
l.U(U.y/ — 1.1 J 


l.U(U.9b— 1.1 j 


l.lU^l.Uo— l.Lo) 





13(0.079-0.16) 


0.081(0.038-0.11) 


1.7(1.1-2.1) 


1.7(1.0-2.5) 


Abundances (solar) 


















O 


d 


d 


d 





21(0.086-0.37) 


0.24(0.058-0.52) 


d 


d 


Nc 


1.0(0.58-1.8) 


0.99(0.45-1.8) 


1.5(0.98-2.2) 





75(0.53-1.1) 


0.48(0.22-0.96) 


0.081(0.015-0.22) 


0.22(0.056-2.1) 


Na 


d 


d 


d 


3 


7(1.4-8.2) 


2.6(<7.6) 


d 


d 


Mg 


0.13(<0.39) 


0.39(0.11-0.71) 


0.23(0.050-0.45) 




d 


d 


d 


d 


Si 


0.18(0.078-0.35) 


0.12(<0.26) 


0.34(0.23-0.50) 




d 


d 


d 


d 


S 


0.37(0.22-0.65) 


0.24(0.071-0.42) 


0.32(0.19-0.50) 




d 


d 


d 


d 


Ar 


0.13(<0.43) 


0.69(0.27-1.3) 


0.54(0.27-0.88) 




d 


d 


d 


d 


Ca 


0.17(<0.55) 


0.24(<0.73) 


0.41(0.10-0.74) 




d 


d 


d 


d 


Fc 


0.11(0.057-0.17) 


0.10(0.051-0.18) 


0.14(0.097-0.22) 





092(0.056-0.16) 


0.12(0.067-0.20) 


0.019(<0.064) 


0.049(<0.18) 


Others 


0.13(<1.1) 


0(<0.93) 


0.13(<0.84) 





14(0.065-0.22) 


0.080(<0.21) 


0.0080(<0.049) 


0.11(0.022-0.34) 


L x (10 30 ergs s _1 )° 


28 


34 


38 


1 


5 


1.9 


22 


30 


rcduced-x 2 (d.o./.) 


0.953(175) 


0.910(150) 


1.03(275) 


1 


12(96) 


1.07(88) 


1.11(67) 


0.931(85) 



a Thc 1-T model for DoAr 21 and the 2-T models for ROXs 21 and ROXs 31 are used. Parentheses indicate the 90% confidence limits. 
b Scc Figure 1. 

c Thc boundary between phases Q and F2 is not clear because of the existence of a rising part (Figure 1), however there arc no significant differences in parameters 
when the boundary is changed. 

d Thcsc arc assumed to be the same as "Other" elements. 
e Absorption-corrected luminosity in 0.5—9.0 keV. 



